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We investigate functions that are exact solutions to chaotic dynamical 
systems. A generalization of these functions can produce truly random num- 
bers. For the first time, we present solutions to random maps. This allows 
us to check, analytically, some recent results about the complexity of random 
dynamical systems. We confirm the result that a negative Lyapunov exponent 
does not imply predictability in random systems. We test the effectiveness 
of forecasting methods in distinguishing between chaotic and random time 
series. Using the explicit random functions, we can give explicit analytical 
formulas for the output signal in some systems with stochastic resonance. We 
study the influence of chaos on the stochastic resonance. We show, theoreti- 
cally, the existence of a new type of solitonic stochastic resonance, where the 
shape of the kink is crucial. Using our models we can predict specific patterns 
in the output signal of stochastic resonance systems. 
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Recently, many outstanding papers [|T}-j^] have stated the importance of hav- 
ing true random models. The best existing pseudorandom number generators 
can yield incorrect results due to the "unavoidable" correlations that appear be- 
tween the generated values On the other hand, there is a great interest in 
random dynamical systems In recent years, there has been much discus- 
sion about the transition to chaos and the way to characterize predictability and 
complexity in these systems There is also strong controversy about the 
existing methods to distinguish chaotic and completely random systems fl8HIO 



In the present paper we investigate explicit functions that are exact solutions to 
nonlinear chaotic maps. A generalization of these functions can produce truly 
random sequences. Even if the initial conditions are known exactly, the next 
values are in principle unpredictable from the previous values. These functions 
cannot be expressed as a map of type X n+ \ = g (X n , X n _i, ...,X n - r +x). Using some 
of these functions we can exactly solve random maps as the following: 

X n +i = / (X n , I n ) , (1) 



where I n is a random variable. We can confirm the result [§,0] that a negative 
Lyapunov exponent does not imply predictability in random systems. We show 
that the forecasting methods []8| — pTOfl are very effective in distinguishing chaos 
from random time series. We investigate the influence of the level of chaos on 
the stochastic resonance (SR). We can give explicit analytical formulas for the 
output signal of some systems with stochastic resonance. We show the existence 
of a new type of solitonic stochastic resonance (SSR), where the shape of the 
kink is crucial. 



I. INTRODUCTION 



There is the common belief that, as truly random numbers should be unpredictable in 
advance, they must be produced by random physical processes such as radioactive decay, 
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thermal noise in electronic devices, cosmic ray arrival time, etc. [|]]. 

Knowing the past and present values should give no information as to future outcomes 
of a truly random variable . Thus, a recursive mathematical algorithm should not be able 
to describe a truly random process. From this, it seems, that deterministic randomness is 
inherently unattainable 0. 

Here we have two problems as a motivation for our work: 

1. How to describe theoretically these physical phenomena that are truly random. 

2. How to produce truly random numbers, which are necessary in different physical cal- 
culations such as Monte Carlo method. 

The purpose of our study is to find explicit functions that produce truly random dynam- 
ics. These functions can be used as random number generators and as analytical solutions 
to nonlinear random systems. 

It is well known |fTI] , |l2]| that the function X n = sin 2 {6n2 n ) is the general solution to the 
logistic map X n+1 = 4X n (1 — X n ). Recently, other chaotic maps have been reported to have 
exact solutions JT3-irS[. In the present paper we will investigate in detail a generalization of 



the solution to the logistic map: 

X n = sin 2 (Ottz 71 ) , (2) 

where z is a real number. 

For z integer, function ([|) is the general solution to the family of maps: 

X n+ i = sin 2 (^z arcsin \j X,^j . (3) 

Even for a real z we can calculate the Lyapunov exponent of map (||) exactly: A = In z. 

For z > 1, map (0) is chaotic. Nevertheless, for fractionary z the dynamics contained 
in function (0) is quite different from that of map (j^). In fact, for a fractionary z, the 
first-return map generated by Eq. (g) is multivalued (see Figs. 1 and 2). Let z be a rational 
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number expressed as z = p/q, where p and q are relative prime numbers. Then the first- 
return map produced by function (Q) is a curve such that, in general, for a value of X n we 
will have q values of X n+ \. On the other hand, for a value of X n+ \ we will have p values 



of X n . Geometrically, these curves are Lissajous figures |T7|]. But we should note that their 
meaning here is very different from that in their original definition. In this context, they 
represent chaotic first-return maps. For z irrational, the first-return map is a random set of 
points as shown in Fig. 3. 

The paper is organized as follows. In Sec. II we study the properties of the functions 
X n = sin 2 (9irz n ). We present a rigorous proof that, for z fractionary, the produced se- 
quences are absolutely unpredictable in advance. Moreover, the outcomes are completely 
independent. In Sec. Ill we discuss the use of these functions in actual numerical calcu- 
lations. Section IV is dedicated to random maps of type X„ +1 = f(X n ,I n ), where I n is a 
random variable. Function (0) can help one to find analytical solutions to these maps. We 
calculate exactly the complexity of a random map. This allows us to check some recent 
results about the complexity and predictability of random maps. In Sec. V we address the 
problem of distinguishing chaos from random time series. For this, we check the effectivity 
of the so-called "nonlinear forecasting methods." Section VI is devoted to stochastic reso- 
nance (SR). First we give some introductory remarks about the historical developments in 
SR. Considering the fact that we can calculate exactly the Lyapunov exponent of a class 
of chaotic maps, we are able to investigate the influence of the level of chaos on SR. This 
is done first in the most common setup for SR: a bistable system. Then, we investigate 
the so-called nonlinear static systems with SR. For these systems, we can present explicit 
analytical functions that describe the output of the system. Using the functions we can 
investigate the actual dynamics of the system. Finally, based on theoretical investigations, 
we show the existence of a new type of solitonic stochastic resonance, where the shape of 
the kink is crucial. 
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II. EXPLICIT STOCHASTIC FUNCTIONS 



After a rigorous analysis of function (Q) we arrive at interesting conclusions. For most 
fractionary z > 1 function (^|) is not only chaotic, but its next value is impossible to predict 
(from the previous values) unless 9 is exactly known. When z is an integer, the initial 
condition A defines univocally the value of 9 (any value of 9 out of the interval < 9 < 1 
defining Xq is equivalent to one in that interval). If z is fractionary, this is not so. There 
exists an infinite number of values of 9 that satisfy the initial conditions. The time series 
produced for different values of 9 satisfying the initial conditions is different in most cases. 
The fact that we know the initial conditions does not imply that we can determine 9. So 
the next value is unpredictable. 

Let us consider the case z = 3/2 (see Fig. 2). If we wish to calculate X n+ i from the 
value A„ we will have two choices: 



X, 



1 



n+l 



l±(l-4X n )(l-X n ) 1 / 2 . (4) 



The value X n+ i could be expressed as a well-defined function of the previous values if 

1 /2 

(1 — X n ) could be a rational function of the previous values. However, each time we try to 
do this we meet the same difficulty because the previous values are also irrational functions 
of the past values. This process can continue up to infinity. 

A different way to see this phenomenon is the following. Consider the family of functions 

^ = sin 2 [(# + £;)7r^], (5) 

where 9 = 9 + k, k is integer. 

For all k, the time series A^ {k fixed, n as time) have the same initial conditions. If 
z is an integer, the initial condition defines the complete sequence (see Table I). However, 
for z fractionary all the time series are different. This is because the period of function A^ 
(now n is fixed and k is variable) is different for different n (for instance, when z = 3/2, 
the period of X% is 2 n ). In general, for z = p/q, the period is q n . That is, A n+1 cannot 
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be determined by X n . Moreover, X n+ i cannot be determined by any number of previous 
values. Let us see the following example with z = 3/2. Suppose X n = 0. Now we have two 
possibilities X n+ \ = or X n+ \ = 1 (see Table II). Assume 9q = and n — 0. For any 9 = k 
(k integer), X n = 0. Now, X n+ \ = sin 2 [(3/2) kit). So, X n+ \ = for k even, and X n+ i = 1 
for k odd. But there is no way we can know k from the statement X n = (for all k integers 
this statement is true). This uncertainty about the next value is present for all points X n 
except X n = 1/4 and X n = 1. But these two points are a set of zero measure. That is, for 
almost all the points in the interval < X n < 1, the next value is unpredictable. 

For z irrational there are infinite possibilities for X n+ \. All values are unpredictable. 
But let us continue with the simple case z = 3/2. Suppose now that 9 = 2 m , where m is an 
integer. Note that in this case X = 0. But, unless we know 9, we never will know when the 
value X m+ i will be equal to 1 (see Table II). We can have a string of m + 1 zeros (m can 
be as large as we wish) and only in the point X m+ i does the sequence change from a string 
of zeros to the value 1. So, for any finite number m + 1 of previous values X , Xi, X 2 ,..., 
X m ; the next value is not defined by the previous values. Note that in this example we can 
have a string of zeros, but this is because the value X n = is a pseudofixed point of the 
map (X n , X n+ i) due to the intersection of the graph in Fig. 2(a) with the line X n+ i = X n . 
However, in general, the sequence is very stochastic. On the other hand, the uncertainty 
about which is the next value remains for all the points in the interval < X n < 1 except 
for X n = 1/4 and X n = 1. The general uncertainty increases for p > q > 2 (see Table III). 
In this case, the unpredictability is true for all values of X n . 

On the other hand, if z is irrational, then the points on the first-return map (X n ,X n+1 ) 
will fill the square < X n < 1; < X n+ i < 1 (see Fig. 3 and Table IV). For a large but 
finite number n, the map is an erratic set of points (we should exclude the numbers of type 
z = m 1 '*, where m and k are integers, because in this case the sequence is predictable given 
k previous values). 

Note that we can consider X^ defined by Eq. ([5|) as an infinite matrix, where the 
"columns" are the stochastic sequences (dependence on n) and the horizontal "rows" are 
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periodic (or quasiperiodic for irrational z) sequences that represent the dependence on k. 
For z = p/q, the "rows" are periodic sequences with period q n (see Tables II and III). We see 
that all the row sequences have different periods. So, all the column sequences are generally 
different. However, for each integer m, there is an infinite set of columns having a string of 
values of length m that is identical in each number of this set. That is, in the matrix 
given an initial string of length m = 2, we will find a string identical to it with a period q 2 . 
Note (in Table II) that the string (0, 1, 1/2) can be found in infinite columns. However, the 
next value is always uncertain. It can be X3 = 0.1464... or X3 = 0.8536... . Just to know that 
the previous values are (0, 1, 1/2) does not give us the knowledge to determine the next value. 
The string (0, 1, 0.5, 0.8535..., 0.0380..., 0.9157..., 0.8865..., 0.0711..., 0.8447..., 0.9686...) can 
be found with period 2 9 . That is, the column number 2 9 + 1 possesses this same string. 
However, the value X w is not always X w = 0.7544... . It can be X w = 0.2455... with the 
same probability. 

In general, given an initial string of length m, we will find a string identical to it with 
a period q m . At the same time, most of these strings possess different next values (we have 
seen a striking example in the above-given text). Suppose there is a univalent function 
X n+ i = g (X n , X n -i, X n - r+ i) that is equivalent to the sequence (fj) for z fractionary. 
If we have more than one sequence Xq, Xi, X2,..., X m _i with different next values, then 
we should decide that the map we are looking for cannot be of order m. If for any m, 
m = 1,2, 3, 00; we have more than one sequence X , Xi, X 2 ,..., X m _i, such that the next 
values are different, then such a map does not exist. 

In the above-given text we have shown that for each string of values X , Xi, X 2 ,..., X m -i, 
there is another sequence with these same values but with different proceeding values. 

For z irrational, all the row sequences are quasiperiodic and different. The column 
sequences correspond to completely random sequences. These functions can produce a set 
of completely independent values. 
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III. RANDOM NUMBER GENERATORS 



Now we should say some words about the use of these functions in actual numerical 
calculations. The argument of function (12) increases exponentially. So, there can be some 
problems in generating very large sequences. A practical solution is to change parameters 
9 after a fixed number n = N of sequence values X n . Suppose N is a number for which 
there are not calculation problems. For producing the new set of values of X n (with a new 
9) we start again with n = 0. This procedure can be repeated the desired number of times 
(remember that even if the sequence is finite, it will be unpredictable; and a sequence formed 
as a set of unpredictable sequences will be also unpredictable). It can be shown that there 
exists always a 9 such that, with it, the original function will produce the same sequence as 
that generated with the procedure of changing 9. 

For the calculation of truly random numbers with function (§) the best way is to use an 
irrational z. This irrational z does not have to be a large number. For instance, we can 
use z = it. The geometrical place of the return map for z irrational is the whole square 
< X n < 1; < X n+ \ < 1. So, we do not have to worry about the method for determining 
the next value of 9. For example, we can use the following method in order to change 
parameter 9 after each set of N sequence values. 

Let us define 9 S = AW S , where s is the order number of 9 in such a way that s = 1 
corresponds to the 9 used for the first set of N values X n ; s = 2 for the second set, etc.; W s is 
a "stochastic" sequence. For instance, the values W s can be obtained from the same sequence 
X n . The inequality A > 1 should hold in order to keep the absolute unpredictability. 

Another important question about good random numbers is to have a generator 
able to produce uniformly distributed points. By means of the transformation Y n = 
(2 /it) arcsin ^X^ 2 ^; we can obtain random numbers uniformly distributed on the inter- 
val (0, 1) ||19|| . Once we have uniformly distributed random numbers, we can use well-known 



transformations to generate random numbers with any given distribution ||T9 



We have performed several standard statistical tests with the functions X n = sin 2 (9iiz n ) 
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(after the transformation Y n = (2 /it) arcsin (^X^ 2 ^j. Among them are the following: the 
central limit theorem test, the moments calculations, the variance calculation, and the \ 2 
test. The sequence Y n has passed all these tests satisfactorily. For instance, the theoretical 
values for the moments and variances are the following: < X n >= -4r> <^ n = ,„ , , "v ,^ 2 , 
and these values are obtained when we use the sequence Y n . 

The autocorrelation function C m =< YiY i+m > — < Yi > 2 (where <> is the average 
overall % with % = 1, 2, 3, ...) can be shown to be zero even for m = 1. For the known chaotic 
maps (which sometimes are used as pseudorandom number generators) \C m \ decays with m, 
but there is a range of this dependence that is related to the correlation or memory time. 

Recently [^TJ a new method has been developed, which allows us to compare the ran- 
domness of different sequences. In these works, a measure of randomness (we will call it R) 
is introduced. 

Suppose we have a sequence of values Ui, U2, U3, U n . Form a sequence of vectors 

= [Ui, U i+1 , U i+m -i] . (6) 



Now, we will define some variables: 



number of j such that d 
C™( r ) = 



N-m + 1 



< r 

(7) 



where d 



is the distance between two vectors, which is defined as follows: 



max ( I t/j+fc-i -Uj+k-i\) (8) 



with k = 1, 2, m. 

Now we can define the measure of randomness: 

R(m, r, N) = <^ - (fffi 1 , (9) 

where 

1 N-m+1 

^ ln^ m (r). (10) 



(r) N-m + 1 ^ 



This measure depends on the resolution parameter r and an "embedding" parameter m. 



This technique has been proved to be very effective in determining system randomness jzT 

For given r and m we have a maximum possible randomness. A sequence with maximum 
randomness is uncorrelated. The randomness of our sequences Y n with z irrational is the 
maximum possible for the given r and m. For instance, if r = 0.025, the maximum possible 
randomness is R = In 40. The randomness of function (g) with z = ir approaches the value 
R = 3.688 for increasing N. For comparison, the randomness of the logistic map at the 
point of full chaos is R = 0.693. Even if we further decrease r and increase m and N, for 
the logistic map and other usually chaotic maps, R saturates and remains constant. 

On the other hand, for r — > 0, the randomness of function (0) with z = it tends to the 
maximum possible value, i.e., R — > ln(l/r). For r — > 0, it never saturates. 

We should say that the pseudorandom number generators described in Ref. 1 can pass 
some of the statistical tests devised to check pseudorandomness f|]. However, hidden er- 
rors in these generators have been found |JJ. Several researchers have traced the errors 
to the dependence in the pseudorandom numbers. Indeed, they are all based on recursive 
algorithms. 

Recently simulations of different physical systems have become the "new tests" for pseu- 
dorandom number generators. Among these systems are the following: the two-dimensional 
Ising model, ballistic deposition, and random walks. Nogues et al || have found that using 
common pseudorandom number generators, the produced random walks present symme- 
tries, meaning that the generated numbers are not independent. On the other hand, the 
logarithmic plot of the mean distance versus the number of steps N is not a straight line 
after N > 10 5 (in fact, it is a rapidly decaying function). 

D'Souza et al. use ballistic deposition to test the randomness of pseudorandom number 
generators. They found correlations in the pseudorandom numbers and strong coupling 
between the model and the generator (even generators that pass extensive statistical tests). 

In a ballistic deposition model of growth, free particles initiated at random positions 
above a one-dimensional substrate, descend ballistically and stick upon first touching the 
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surface of the growing cluster. The substrate of length L consists of discrete columns indexed 
by integer values of x with 1 < x < L. The growth interface is defined by the maximum 
occupied site along each column h(x,L), where h(x,L) also takes on discrete values. 

The width of the growth interface on average increases following a power law 

behavior until reaching a steady asymptotic value, the magnitude of which depends on the 
underlying substrate size L. For £l(£) we have 

e L (t) = jj:[h(x,t)-<h(t)>]\ (ii) 

x—l 

where < h(t) > is the mean height of the surface at time t. 

One consequence of the Kardar-Parisi-Zhang theory is that the steady state behavior 
for the interface fluctuations in one dimension should resemble a random walk, i.e., — > 
oo) ~ L 1 ! 2 . Thus, a random walk again serves as a good test for random numbers. These 
two papers pm have been cited by Fisher || in his article about the great problems of 
statistical physics for this century. 

With our numbers Y n , the produced random walks possess the correct properties, includ- 
ing the mean distance behavior < d 2 >~ iV (see Fig. 4). 



IV. RANDOM MAPS 

The functions of (Q) are important not only for numerical simulations. They are relevant 
by themselves as theoretical paradigms of stochastic processes ||. Considering the fact 
that these are explicit functions, we can use them to solve (analytically) many theoretical 
problems in stochastic dynamical systems. 

Consider the following random map ||: 

1 



X 



n+l 



l + 4(l-4X n )(l-X n ) 1/2 , (12) 



where I n is a random variable that takes the values ±1 with equal probability. 
An exact solution to this random map can be written as follows: 



11 



X n = sin 2 [9n (3/2) n ] . (13) 

Now we can check some of the results discussed by the authors of Refs. 6 and 7. They 
have introduced a measure of complexity K in terms of the average number of bits per time 
unit necessary to specify the sequence generated by the system. In dynamical systems of 
type ([!]) [Eq. (|1^ is an example] this measure coincides with the rate K of divergence of 
nearby trajectories evolving under two different realizations of the random variable I n . 

The complexity of the dynamics can be measured as 

K = \6 (A) + h, (14) 

where A is the Lyapunov exponent of the map, h is the complexity of I n , and 6(X) is the 
Heaviside step function. Complexity h should be defined also as the average number of bits 
per time unit necessary to specify the random variable I n . When I n is a usual chaotic noise, 
then h coincides with the Kolmogorov-Sinai entropy. 

In the case of the random map (|T2j) A = In (3/2) and h = In 2. Hence, K = In 3. On the 
other hand, any calculation (theoretical or numerical) of K for the dynamics generated by 
function (|T^) yields the correct value K = In 3. Moreover, even an independent calculation 
of the complexity of this dynamics using different methods pO^pif produces the same result. 

Using the functions (0) we can also solve the map 

X n+1 = 1 + (15) 

where I n is defined as in Eq. (|12"D. 

Here the Lyapunov exponent is negative: A = In (1/2) < 0. However, the complexity is 
positive: K = In 2. 

In the presence of random perturbations, K can be very different from the standard 
Lyapunov exponent and, hence, from the Kolmogorov entropy computed with the same 
realization of the randomness. 

We stress that a negative value of A does not imply predictability. 
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In general, if we apply the measure of complexity K to our function (0), then we obtain 
the following results: For z = p/q (p and q relative primes), K = Inp. If z is irrational, the 
complexity is infinite!. 

We should say that using function (Q) we can create complete sets of orthogonal elements. 
In the same way that we can solve, to begin with, any map of type X n+ i = f (X n ), we can 
also solve theoretically many important problems in stochastic dynamical systems. 



V. NONLINEAR FORECASTING METHODS 

Now we address the problem of deciding which of the proposed methods p HlQfl for dis- 
tinguishing chaos from random time series are more effective. Recently a new method based 



on nonlinear forecasting was proposed |§-10|. The idea of the method is as follows. One 



can make short-term predictions that are based on a library of past patterns in a time series 
(the method of nonlinear forecasting is described in Refs. 8-10 and the references quoted 
therein). By comparing the predicted and actual values, one can make distinctions between 
random sequences and deterministic chaos. 

For chaotic (but correlated) time series, the accuracy of the nonlinear forecast falls off 
with increasing prediction-time interval. On the other hand, for truly random sequences, 
the forecasting accuracy is independent of the prediction interval. The decrease with time 
of the correlation coefficient between predicted and actual values can be used to calculate 
the largest positive Lyapunov exponent of the time series §| . 

Function (|2|) is a very good model system to check this and other methods. In fact, 
for z integer, these are chaotic sequences of type X n+1 = / (X n ). For z = m 1//fc , we have 
chaotic maps of type: X n+ i = g (X n ,X n -i, ...,X n -k+i). For z fractionary, we have different 
types of random sequences with different complexities. Finally, for z irrational (generic), 
the sequence is maximally random. 

Suppose we have a sequence U\, U2, Un- Now we construct a map with the dependence 

jjpredicted ag & "f unction » Q f U observed _ 
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If we have a correlated chaotic sequence, this dependence is a straight line, i.e., 
jjpredicted ^ jjobserved ^jjgjj ^he forecasting method is applied for one time step into the 
future). When we increase the number of time steps into the future, this relation dete- 
riorates. When we apply this method to function (|2]) with z = n [after transformation 
Y n = (2/tt) arcsin \J X n ] , even the map jj^ edlcted vs U° bserved for one time step into the future 
is a map equivalent to that shown in Fig. 5. On the other hand, the correlation coefficient 
is independent of the prediction time. In fact, there are no correlations. The details will be 
given elsewhere. Nevertheless, we should say that this method is quite efficient in distin- 
guishing chaos from randomness. However, it cannot distinguish between different random 
time-series. 

Other methods discussed in Ref. 8, are less effective in this task. They are more qualita- 
tive, requiring subjective judgment about whether there is an attractor of given dimensions. 



VI. STOCHASTIC RESONANCE 



A. Lyapunov exponents 



A phenomenon that has awakened very much interest in last several years is stochastic 
resonance (SR) [f2"2] -4Tj|. The classical model for this phenomenon is the following: 



x — x + x 3 = Aq sin (cut) + 77 (t) . (16) 

The sum of a noise signal, rj (t), and a weak periodic signal is used to drive a bistable 
system. The most important characteristic of SR is that the signal-to-noise ratio (SNR) has 
a maximum in the plot SNR vs D, where D is the noise intensity, for a finite nonzero value 
of the noise intensity. 

It has been shown that SR still occurs when chaos, rather than noise, is used as the 



nonperiodic component of the driving signal f2~7fl . Several authors have investigated the SR 



in chaotic systems 27-32 . 
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A very important question is how SR depends on the largest Lyapunov exponent of the 
driving chaotic noise p5] . We address this problem systematically for the first time, since 
we can solve exactly the problem of calculating the Lyapunov exponent. We should say that 
for large values of the Lyapunov exponent, the SR is not a very sensitive phenomenon on 
the level of chaos (this is unlike the phenomena discussed in Refs. 1-3. The curve SNR vs 
D practically has no variation for A 3> In 3. However, in the interval In (3/2) < A < In 2, 
the SR strongly depends on the Lyapunov exponent. In this interval, the maximum of 
SNR is shifted to the right (larger noise intensities) and is amplified!. Figure 6 shows 
the dependence SNR(D) for different values of the Lyapunov exponent A. These data are 
the result of numerical simulations of Eq. (0). We can compare this result with that 
obtained in Ref. 33. In this work, the phenomenon of stochastic resonance is studied in the 
presence of colored noise. In overdamped systems, the authors find that SR is suppressed 
with increasing noise color. In contrast, for colored noise induced by inertia (as well as for 
asymmetric dichotomic noise), they obtain an enhancement of SR. 

The same result can be obtained in the so-called threshold systems |3~3|-|3"6] . For instance, 
define 

In = g(Pn + Vn), (17) 

where p n is a periodic function, r] n is some kind of noise, and g (x) is a function with some 



properties that allow the existence of SR I34H36 
The simplest case is the following: 

( -V,x < x th , 

g(x) = \ (18) 
[V,x> x th - 

A nonlinear circuit with this kind of threshold nonlinearity is discussed in Ref. 36. 

Different measures have been used to characterize stochastic resonance. In particular, in 
Ref. 41 the dynamics of noisy bistable systems is analyzed by means of Lyapunov exponents 
and measures of complexity. 

It can be shown that, in stochastic resonance systems, the function K (D) (where K is 
the complexity as defined previously and D is the noise intensity) has a local minimum for 

15 



a nonzero value of D. This minimum represents an optimal value of noise intensity at which 
SR occurs. This result confirms the findings of Ref. 41. 

We can consider K as a "dynamical measure" of SR because K coincides with the rate 
of divergence of nearby trajectories evolving under two different noise realizations. 

This measure can be used to characterize more general stochastic dynamical systems 
such as the following: 

X n+1 = F(X n ,p n ,rj n ), (19) 

where p n is a periodic function and r\ n is the noise. 
Let us investigate a particular example: 

X n+ i = cos [ 1 + e (p n + Drjn) 2 arccos (X n )j . (20) 

Here p n is a periodic function of amplitude a and r\ n is a chaotic noise defined as follows: 

r] n = Y n - 5 - 1, (21) 

where 5 is a parameter for which < 5 < 1 and Y n+ i = sin 2 [z arcsin y/Y^\ . 

For this dynamical system, function K (D) has a minimum for a finite D. 

Let us suppose that p n is a period-one function. We can write down an analytic expression 
for K: 

K = \n(z \ + £ -{{a- D5f + (a- D(z-l)f) }. (22) 

The function K (D) is shown in Fig. 7. For a fixed z, the minimum of K is obtained 
approximately for D = a/ (z — 1). If we minimize K with respect to both D and we 
obtain D = a/5, z = 1 + 5. 

Note that for 1 + 5 < z < 3, the minimum is deeper and is shifted to the right as z is 
decreased. This is a phenomenon similar to that obtained using SNR in Fig. 6. See also 
subsection VI E, where some experiments are mentioned. 

We should say that K can characterize the dynamics of dynamical systems of type ( |T9"D 
even when there is no periodic function at all. 
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In some cases, for a finite value of D we can find the least complex dynamics. 

In some sense, this is a more general phenomenon than the usual stochastic resonance. 
In fact, this is an example of the so-called noise-induced disorder-order transitions, of which 
the SR phenomena can be a subset. 

In Fig. 8, different time series and return maps produced by the dynamical system 
( P0|) are shown. Note that for some intermediate value of D we obtain the least complex 
dynamics. 

In Fig. 9 we see that when both D and z are very close to the optimal values, then the 
resulting dynamics is very predictable. 

Note that this system can be chaotic even when D = 0, due to the intrinsic nonlinear 
dynamics of the system. However, for some finite value of the noise intensity D > 0, we can 
control this chaotic dynamics. So, in this case, we are truly controlling the chaotic system 
using chaotic noise. 



B. Explicit output functions for SR systems 

Recently scientists have learned that stochastic resonance can appear not only in bistable 
systems [[36H3811. A very interesting class of systems is that of the so-called nonlinear static 
(or "nondynamical" ) systems. In Refs. 36-38 a theory of these systems is presented. Using 
this theory and our function (|2|) we can write down an explicit "solution" function to these 
systems. For instance, the function: 

I n = tanh {B [A sin (wn) + D cos {6nz n ) - V th ] } (23) 

can behave as a SR system. Figure 10 shows that the function (^) is a SR system for B = 24. 
For B = 1 the stochastic resonance disappears. Here the SNR was calculated numerically 
using function ( |23"D as the output signal. In fact, we can construct a very general class of 
SR functions of type I n = g (V n ) + £ n , where V n = p n + r\ n is the input and I n is the output. 
Function represents the intrinsic noise [B3 . 
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Note that although the systems described in Refs. 36-38 are called "static systems," 
once we have constructed our explicit functions, e.g., Eq. (p3|), we can obtain exact solu- 
tions to very dynamical systems. Also note that it is very easy to check that there is a 
maximum in the dependence SNR vs D\ however our explicit functions can be investigated 
using mathematical analysis, not only statistics. Very different functions with very differ- 
ent dynamics can have the same SNR behavior and other statistical properties. Using our 
explicit functions we can investigate the true dynamics of the system. 

The analysis of function (fj) allows us to construct a continuous and differentiable function 
with properties similar to those of the chaotic functions Let us give an example: 

f (t) = sin{i?i sinh \a\ cos (coit) + a 2 cos (uint)] + 

(24) 

B 2 cosh [03 cos {uj^t) + a 4 cos (u^t)]}. 
Using functions of this kind we can find analytic solutions to continuous chaotic dynam- 
ical systems. Function (p4|) with the parameter values B\ = 20, B 2 = 30, a\ = 10, a 2 = 
15, as = 10, 04 = 15, lui — 1, u>2 = 7r, cl> 3 = a/2, uj a = e; behaves as a chaotic system (see Fig. 
11). Any investigation (theoretical or numerical) will give the same result: The maximum 
Lyapunov exponent is positive. Moreover, if we need a continuous dynamics with a chaotic 
Gaussian-like "noise" we can use a transformation of Eq. (p4]): g(t) = In [f 2 (t)/ (1 — f 2 (t))]. 
We have been able to produce SR with function (|54T ), f(t), and g{t). In Fig. 12 SNR is 
calculated from numerical simulations of Eq. (|T^) and using the continuous chaotic function 
g(t) = In \p{t)/ (1 - /*(*))], where /(*) is defined by Eq. @. 

C. Solitonic stochastic resonance 

The spatiotemporal SR in the <£> 4 model has been considered in a very interesting paper 
p9| . Recently we introduced the concept of solitonic stochastic resonance (SSR) [0 where 
a soliton moves in a bistable potential created by space-dependent external forces driven by 
a periodic signal and noise. This seems to be equivalent to the conventional setup for SR, 
however the conditions for the existence of SSR are different from that of SR with a point 
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particle in a bistable potential. The situation for SSR can produce very interesting phenom- 
ena like the transformation of the soliton into a three- "particle" system of two solitons and 



an antisoliton IHO 



Here we will present another framework for SSR. 
The function 



(x, t) = tanh {B [x - x - A sin (cut) - Df (t)} } , (25) 



where f{t) is defined in Eq. fl24|), can be used to find analytic solutions to nonlinear partial 
differential equations. 

Recall that function ( f25|) is a SR solution [see Eq. fl23")1. For instance, if we take the 



time series produced by (ft {x = 0,t) [with B = 12, A = 0.67, to = 0.88, x = 2, and f(t) is 
the function (|23|)] we obtain a new kind of SSR from (|25|). In fact, using the solution ( |25| ) 
we can prove that the overdamped perturbed ip 4 equation 



B 



uAcos (ujt) + Df(t) 



cosh [B (x — xo) 

where F(x) = a tanh [B (x — x )), possesses a different kind of SSR. 

We can calculate analytically the SNR for the dynamics of Eq. (|26|) . 

Suppose that, in Eq. (|26|), instead of f(t), the noise is described by the function h(t) = 
(2/7r) arcsin/(t), and f(t) is given by function (|24|). This is equivalent to a uniformly 
distributed noise in the interval (0, 1). We will define b/2 = xo — A. 

Following the ideas of Ref. 34 we can calculate an approximate analytical expression for 
the SNR when B > 1: 

' for f < | - A 
SNR =1 i(f - | + 2A) for § - A < f < § + A 
& far^>i+A 
It is evident that there is a maximum in the curve SNR(D). 
Our theoretical results on the theory of solitons perturbed by external forces 



allow us to understand the dynamics of Eq. ( P6|) and to interpret the physics of solution 
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(f2~5|). In order to obtain the desired dynamics we should solve Eq. fl2~6|) with an initial 
condition representing a soliton situated in a vicinity of point In this case, the 

soliton center of mass will be oscillating inside the potential well created by the force F(x). 
This is exactly what represents solution (|25|). 

When we investigate the time series obtained after the numerical simulation of Eq. 



we obtain SSR as predicted by the theoretical solution (j25|) . 

The SNR vs D plot depends on the value of B. For very large values of B, the SNR 
vs D plot has a very nice maximum (see Fig. 13). The SNR(D) dependence shown in Fig. 
13 was calculated numerically from the time series generated by the function <p(x = 0, t) as 
a solution of Eq. (^). The same result is obtained if we investigate the analytic solution 
(p5p. For very small values of B, the SSR disappears. Thus, this is a SSR that depends on 
the shape of the soliton. In particular, it depends on the width of the soliton, which can be 
expressed as S = 1/B. This SSR is different from the one obtained for a soliton moving in 
a bistable potential well [[40] and that described in Ref. 39. 

In Ref. 39 the synchronization of a linearly coupled chain of N overdamped bistable 
elements, subject to a deterministic periodic signal and uncorrelated white noise, is addressed 
in the continuous limit of a ip A theory. The cooperation between noise and coupling is shown 
to lead to spatiotemporal stochastic resonance. There, the bistability of the ip 4 equation 
on the potential U((p) ~ (if 2 — l) 2 plays the most important role. On the other hand, in 
our previous paper j|(J we considered the stochastic resonance of a soliton moving in a 
bistable potential created by inhomogeneous external forces F(x). In this paper, the output 
signal is the coordinate of the soliton center of mass. In the present work, the relevant 
output signal is <p(x = 0,t). The soliton is moving in a monostable well potential created by 
inhomogeneous external forces. However, the most striking feature is that the width of the 
soliton determines the existence or not of the solitonic stochastic resonance. The shape of 
the output signal possesses patterns that are very different from that obtained in a bistable 
system. They are more similar to the patterns that appear in threshold systems. 
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D. Patterns in the output signal of SR systems 



Once we have an explicit solution that describes the stochastic resonance system as the 
following: 



where P(t) is a periodic function, r)(t) and £(t) are different manifestations of noise dynamics, 
we can calculate the SNR exactly. We should note that SNR is the main measure of stochastic 
resonance and is widely used in SR literature. 

For instance, let us define the different noises as follows. is a white noise with zero 
mean and correlation function: 



where Q is a constant parameter. r)(t) is a Gaussian noise with zero mean and correlation 
function: 



I(t)=g(P(t)+ V (t)) 



(27) 



and 




(28) 




(29) 




(30) 



Following Ref. 38 we can consider the case 



g(v) = v\ 



(31) 



where 




(32) 



In this case the SNR is 



SNR = ir 



18aV 4 + 9aV 2 + (9/8)a 6 



(33) 



2Q + t f (44a 6 + 54aV 4 + (27/2)aV 2 ) ' 
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This is exactly the SNR obtained in Ref. 38. 

The static character of the present nonlinearities allows a direct statistical analysis, in 
which all quantities relevant to characterize the SNR in the output signal can be obtained 
from statistics computed directly on the input noises. In fact the SNR is a statistical measure 
based on the statistical properties of noises r)(t) and £(£). 

Nevertheless, we believe that using our explicit functions we can obtain much more 
information about the output signal. Some of this information can have statistical character, 
but we will have also dynamical and geometrical information about the output signal. 

For instance, we can predict the values of the local maxima and minima in the time 
series, and the distance between them. We can obtain the exact analytical shape of the 
extrema. 

In any stochastic resonance output signal there are patterns. These patterns can be 
different for different systems. The following function 



tanh[B(Asm(wn) + D cos(0irz n )) - V th ] + 1 

2 



(34) 



is the analytical solution for a circuit with I — V characteristic of type: 




(see Ref. 36). A typical time series is shown in Fig. 14(a). 
In Fig. 14(b) is shown the typical time series for 




-1 for V < -0.5 



for - 0.5 < V < 0.5 



(35) 



1 for V > 0.5 



In Fig. 14(c) we show the function 




(36) 



where 




(37) 
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which is the output signal for a circuit system with the characteristic 

( for V < V th 

g(V) = \ . (38) 

{ b(V - V th ) for V > V th 

Figure 14(d) shows the output signal 

i n = 9(V n ) + { n , (39) 

where g(V) = V 3 , V n = asm(wn) + rj n , rj n = DY n , Y n = ln(j^), X n = sin 2 (6-Kz n ), 
£ n = Qlrij^-, and Z n = sin 2 (9nn n ). 

All these systems present stochastic resonance. All these systems can be tuned to have 
the same SNR. However, note that all the patterns are different. Compare them to the 
typical time series of a classic bistable stochastic resonance system shown in Fig. 15. The 
information about these patterns is in the explicit functions that can be written down using 
our stochastic functions. 

We can even make predictions about the outcomes in these stochastic systems. For 
instance, in the function ( p4|) we can say that with a probability p = 0.8, after the function 
I n has taken the value I n — 1, it will take the value I n = 0. Meanwhile, we can expect that 
it will remain in the state I n = for an average time aproximately equal to the period of 
the periodic input signal. 

On the other hand, the output function (|36| ) will give us much more information about 



the actual shape of the input periodic signal than the functions (|34|) and (p5|) . 

Systems with intrinsic and external noises are expected to be very random. Nevertheless, 
using the theoretical information obtained from our explicit function (|39|), we can make very 
remarkable predictions. For instance, if the output signal takes a "large" negative value (say 
I n = —20), then with absolute certainty we can predict that the next values will be negative 
and \I n \ — ► 0. When I n reaches the value I n = 0, we can predict that the next value will be 
positive, but the exact value is unpredictable. When it takes this positive value, the next 
value will be negative with absolute certainty. The larger the absolute value of I n when 
it takes a positive value, the larger the absolute value of the next negative value that it 
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will take. Note that all the randomness of this dynamics is produced when I n is near zero. 
When I n is far from zero, we can make exact predictions of the next values. All this can be 
corroborated when we observe the first-return map of this dynamics (Fig. 16). 

We can see the stochastic resonance as a phenomenon that transforms a complex dynam- 
ics into a simpler one. That is, the output signal is less complex that the input signal. But 
there are also phenomena that lead to a more complex behavior (e.g., the chaotic systems). 

Using our functions we can predict the existence of new complex phenomena. 

After an analysis of the functions X n = sin 2 (9irz n ), which we have shown to produce 
complex dynamics, the first characterics that surface are the following: The function can be 
rewritten in the form X n = h(f(n)), where the argument function f(n) grows exponentially 
and the function h(y) is always finite and periodic. 

However, a more thorough analysis shows that (to produce complex behavior) the func- 
tion f(n) does not have to be exponential all the time, and the function h(y) does not have 
to be periodic. 

In fact, it is sufficient that the function f(n) be a nonperiodic oscillating function where 
there are repeating intervals with finite exponetial behavior. For instance, this can be a 
chaotic function. On the other hand, function h(y) should be noninvertible. In other words, 
it should have different maxima and minima. The inverse "function" of h(y) should be 
multivalued. 

The complexity of the output dynamics is proportional to the number of extrema of 
function h(y). 

For example, the following system can produce a dynamics similar to that obtained with 
our function X n = sin 2 (9irz n ): 




(40) 




(41) 



Z n +i — 9{X n ), 



(42) 
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where g(X n ) is a function with several maxima and minima. The first return map is shown 
in Fig. 17. In fact, this is a completely new chaotic phenomenon because the dynamics is 
completely unpredictable. So, when the input is a simple chaotic signal and the system is 
an electronic circuit with the I — V characteristic shown in Fig. 18, then we will have a very 
complex output. This phenomenon is the opposite to the stochastic resonance. Compare 
two I — V characteristic curves for a phenomenon that simplifies the dynamics and for a new 
phenomenon that makes the dynamics extremely complex in Fig. 18. In Ref. 46 a theory 
of nonlinear circuits is presented. There we can find different methods to construct circuits 
with these I — V characteristic curves. 

All the results presented in Sec. VI D, which are related to the Figs. 14-18 were obtained 
through theoretical calculations. 

E. Applications in real systems 

In this section we will present some examples that show how our technique can be used 
in real world applications. 

Our group have designed and constructed a nonlinear circuit (using a concave resistor) 
with the I — V characteristic described by Eq. (36) (see Ref. 46). We wished to check our 
theoretical results about the dependence of SNR on the Lyapunov exponent of the chaotic 
noise. We also desired to observe the patterns for the output signal predicted by our theory. 

In order to have different driving chaotic signals, we produced numerical time series using 
the exactly solvable map ([3]) for different z. 

Then, we transformed the numerical time series into analog signals using a converter. 
These analog signals plus a subthreshold periodic signal were introduced as the voltage to 
the concave resistor circuit. The current in the concave resistor was taken as the output 
signal. 

We should say that the amplification of the SNR in the interval In (3/2) < A < In 2 was 
clearly observed. The maximum SNR for A = In (3/2) was 5 times larger than for A > 3 or 
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for any random noise. 

We will present further details of these experiments elsewhere. 

Thus, in the case that SR is used for the amplification of small signals (as in the dithering 
effect) such that the added external noise is a controllable and manipulable parameter, our 
recommendation is to utilize a uniform chaotic noise with a Lyapunov exponent of the order 
of A « In (3/2). 

Using the techniques described in Ref. 46 it is possible to construct nonlinear circuits 
with the I — V characteristics shown in Fig. 18. 



For our experiments we used the twin-transistor circuit [fig] . 

As an input signal (voltage), we introduced an analog chaotic signal previously produced 
by a nonlinear map. 

Playing with different parameters we were able to produce different unpredictable dy- 
namics very similar to those obtained from our function (Q) and the Figs. 2 and 3. The 
same results can be obtained when we take the input signal from a chaotic electronic circuit. 

In general, the patterns or absence of patterns (it depends on the nonlinearity) predicted 
by our theory are completely confirmed by the experiments. 

Some of these experiments are complicated and should be explained in a separate paper. 

In many relevant applications it is important to get resonances or (in other cases) to 
avoid resonances. 



Recently a new resonance concept was introduced: the geometrical resonance ||7|,fl0 • 
In the usual linear resonance phenomena, the amplitude and frequency of the driving 
force are the most important characteristics. However, in nonlinear systems the shape of 
the driving signal becomes crucial. The geometrical resonance considers the amplitude, the 
frequency, and the shape of the perturbation. 

Suppose we have a nonlinear system A which for some application should be driven by 
a specific driving signal with a given shape. In that case, we can use the output signal of 
another system B as the input signal of the system A. So it is important to predict the shape 
of the output of nonlinear systems. With such information we can design the appropriate 
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system to produce the desired signal needed for the driving of the system A. 

In Sec. VI D we have shown that we can predict specific patterns and the shape of the 
output signal of nonlinear stochastic systems. This is a step forward in the control of chaotic 
and stochastic systems. 

In general, it is very important to find patterns and regularities in the stochastic dy- 
namical systems. In fact, not everything that can be observed can be predicted, only the 



regularities in the observations are the "province of science" |y 



There are many stochastic systems (including systems presenting stochastic resonance) 
where predictions are crucial. In this case we mean predictions of the true values of the out- 
comes using the previous values. Among the concerned areas are the following: geophysics, 



meteorology, climatology, social sciences, etc. [f49|-|52jl 

In Sec. VI D we have investigated a system with external noise and intrinsic noise. These 
noisy perturbations are very unpredictable functions. However, we have shown theoretically 
that we are able to make direct predictions of the time-series outcomes. There are many 
systems with this behavior [38]. When we observe Fig. 14(d), we note that there are very 
remarkable bursts in the time series. In all the mentioned applications it is very important 
to predict these bursts. We have shown that we can do this. 



VII. CONCLUSION 

In conclusion, we can construct functions that are exact solutions to chaotic dynamical 
systems. Moreover, we have generalized functions that cannot be generated by a finite 
recursive algorithm. They can be utilized as theoretical paradigms of stochastic processes. 
Considering the fact that these are explicit functions, we can use them to solve (analytically) 
many theoretical problems in stochastic dynamical systems. Thus, we can apply dynamical 
concepts to describe processes that are usually considered only statistically p8| . 
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TABLE I. Representation of the matrix given by Eq. (5) with z = 2 and 9 = 2 1 / 2 - 1. Note 
that if we start with the same initial conditions, then we will have the same chaotic sequences. 
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0.5151 


0.5151 


0.5151 


0.5151 


0.5151 


0.5151 


0.5151 


0.5151 


x & 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


0.9990 


x 7 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


0.0036 


x 8 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


0.0146 


x 9 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


0.0578 


Xio 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 


0.2181 
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TABLE II. Representation of the matrix X% denned by Eq. (5) with z = 3/2 and 6q = 1. 
Note that all the column sequences possess the same initial conditions Xq = 0. However, all the 
sequences are different in general. 





Oo 


#0 + 1 


6»o + 2 


6<o + 3 


6*0 + 4 


6*0 + 5 


6*0 + 6 


00 + 7 


6*0 + 8 


0o + 9 


00 + 10 


6*o + 11 


6*o + 12 











































Xi 


1 





1 





1 





1 





1 





1 





1 


x 2 


1/2 


1 


1/2 





1/2 


1 


1/2 





1/2 


1 


1/2 





1/2 


x 3 


0.8535 


1/2 


0.1464 


1 


0.1464 


1/2 


0.8535 





0.8535 


1/2 


0.1464 


1 


0.1464 


x 4 


0.0380 


0.1464 


0.3086 


1/2 


0.6913 


0.8535 


0.9619 


1 


0.9619 


0.8535 


0.6913 


1/2 


0.3086 


x 5 


0.9157 


0.3086 


0.4024 


0.8535 


0.0096 


0.9619 


0.2222 


1/2 


0.7777 


0.0380 


0.9903 


0.1464 


0.5975 


x 6 


0.8865 


0.4024 


0.2643 


0.9619 


0.0215 


0.7777 


0.5490 


0.1464 


0.9975 


0.0842 


0.6451 


0.6913 


0.0590 


x 7 


0.0711 


0.2643 


0.5245 


0.7777 


0.9519 


0.9975 


0.9016 


0.6913 


0.4266 


0.1828 


0.0292 


0.0096 


0.1295 


X 8 


0.8447 


0.5245 


0.1213 


0.9975 


0.1923 


0.4266 


0.9087 


0.0096 


0.7674 


0.6214 


0.0649 


0.9784 


0.2751 


x 9 


0.9686 


0.1213 


0.7410 


0.4266 


0.3964 


0.7674 


0.1020 


0.9784 


0.0009 


0.9571 


0.1421 


0.7137 


0.4571 


Xio 


0.7544 


0.7410 


0.0002 


0.7674 


0.7275 


0.0009 


0.7803 


0.7137 


0.0021 


0.7928 


0.6998 


0.0037 


0.8051 



32 



TABLE III. Representation of the matrix defined by Eq. (5) with z = 4/3 and 6* = 1/6. 
Note that the horizontal row sequences possess periods 3 ra . All the next- values in the column 
sequences are unpredictable. 





Oo 


#0 + 1 


6»o + 2 


6<o + 3 


6*0 + 4 


6*0 + 5 


6*0 + 6 


0o + 7 


6*0 + 8 


00 + 9 


00 + 10 


6*o + 11 


6*o + 12 




1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


Xi 


0.4131 


0.9698 


0.1169 


0.4131 


0.9698 


0.1169 


0.4131 


0.9698 


0.1169 


0.4131 


0.9698 


0.1169 


0.4131 


x 2 


0.6434 


0.0531 


0.2014 


0.8431 


0.9177 


0.3019 


0.0134 


0.5290 


0.9966 


0.6434 


0.0531 


0.2014 


0.8431 


x 3 


0.8951 


0.4515 


0.1712 


0.9996 


0.1430 


0.4903 


0.8702 


0.0015 


0.8139 


0.5676 


0.0932 


0.9906 


0.2333 


x 4 


0.9929 


0.6920 


0.2118 


0.0006 


0.2556 


0.7387 


0.9989 


0.7933 


0.3138 


0.0081 


0.1616 


0.6309 


0.9780 


x 5 


0.6475 


0.0675 


0.1584 


0.7792 


0.9668 


0.4302 


0.0018 


0.3463 


0.9291 


0.8462 


0.2261 


0.0308 


0.5633 


X 6 


0.0393 


0.9703 


0.2695 


0.3683 


0.9239 


0.0086 


0.7980 


0.5534 


0.1235 


0.9998 


0.1420 


0.5262 


0.8194 


x 7 


0.4955 


0.5309 


0.4425 


0.5837 


0.3901 


0.6355 


0.3390 


0.6858 


0.2897 


0.7340 


0.2427 


0.7796 


0.1987 


X 8 


0.7550 


0.7849 


0.8132 


0.8400 


0.8651 


0.8884 


0.9097 


0.9290 


0.9461 


0.9609 


0.9735 


0.9837 


0.9914 


x 9 


0.4054 


0.9858 


0.1903 


0.2718 


0.9995 


0.3124 


0.1565 


0.9733 


0.4495 


0.0686 


0.9093 


0.5907 


0.0151 


Xio 


0.0160 


0.6018 


0.9938 


0.4460 


0.0008 


0.5054 


0.9996 


0.5430 


0.0045 


0.4089 


0.9866 


0.6383 


0.0267 
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TABLE IV. Representation of the matrix X% given by Eq. (5) with z = ir and 6*o = 1/4. Note 
that it is difficult even to find "clusters" of equal values in different column sequences. All column 
sequences are completely random and different. 





Oo 


6*0 + 1 


6»o + 2 


6<o + 3 


6<o + 4 


6*0 + 5 


6*0 + 6 


6*0 + 7 


6*0 + 8 


6*0 + 9 


6*o + 10 


6*o + 11 


6*o + 12 


A 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


1/2 


Ai 


0.3897 


0.0516 


0.0457 


0.3761 


0.7983 


0.9995 


0.8307 


0.4169 


0.0646 


0.0348 


0.3494 


0.7756 


0.9976 


A 2 


0.9895 


0.7599 


0.3653 


0.0562 


0.0286 


0.3002 


0.6984 


0.9708 


0.9444 


0.6359 


0.2412 


0.0107 


0.0906 


A 3 


0.4950 


0.4753 


0.4556 


0.4360 


0.4165 


0.3972 


0.3780 


0.3589 


0.3401 


0.3216 


0.3033 


0.2853 


0.2677 


A 4 


0.7996 


0.4643 


0.2603 


0.9388 


0.0012 


0.9002 


0.3252 


0.3937 


0.8535 


0.0114 


0.9684 


0.2003 


0.5356 


A 5 


0.9997 


0.9940 


0.9807 


0.9601 


0.9324 


0.8982 


0.8579 


0.8120 


0.7615 


0.7069 


0.6492 


0.5892 


0.5278 


A 6 


0.7869 


0.5423 


0.1479 


0.9978 


0.0881 


0.6342 


0.7058 


0.0498 


0.9849 


0.2060 


0.4662 


0.8458 


0.0030 


A 7 


0.0521 


0.8342 


0.7685 


0.0216 


0.4881 


0.9847 


0.2517 


0.1484 


0.9368 


0.6171 


0.0003 


0.6509 


0.9186 


A 8 


0.1640 


0.7578 


0.3299 


0.5757 


0.5214 


0.3822 


0.7096 


0.2064 


0.8663 


0.0746 


0.9681 


0.0066 


0.9997 


A 9 


0.5777 


0.8516 


0.9930 


0.9485 


0.7348 


0.4326 


0.1558 


0.0087 


0.0469 


0.2559 


0.5569 


0.8365 


0.9891 


Aio 


0.0013 


0.0343 


0.1084 


0.2171 


0.3507 


0.4976 


0.6446 


0.7789 


0.8885 


0.9639 


0.9982 


0.9885 


0.9357 
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FIG. 1. One- valued first-return map produced by function (|2|) with z = 5. 



FIG. 2. Multivalued first-return maps produced by function (|2|): (a) z = 3/2; (b) z = 8/5. 

FIG. 3. Random first-return maps for z = ir: (a) first-return map produced by function (Q); (b) 
first-return map produced by function @ and with transformation Y n = (2/tt) arcsin 

FIG. 4. Mean distance vs number of steps for a random walk generated with the random numbers 
Y n = {2 /it) arcsin ^X^ 2 j , where X n is given by function (|2|). 

FIG. 5. Predicted values one step into the future vs observed values for the time series generated 
by function (Q) with z = e, after the transformation Y n = (2/tt) arcsin 

FIG. 6. In the interval In (3/2) < A < In 2, the maximum SNR is shifted to the right and is 
amplified. 

FIG. 7. Function K (D) as given by Eq. (22). Solid line, z = 1.5; dotted line, z = 2; dashed 
line, z = 2.5; dot-dashed line, z = 3. Note that as z is decreased, the minimum of function K (D) 
is deeper and is shifted to the right. 

FIG. 8. Time-series and first-return maps generated by the noise-driven dynamical system (20). 
(a) and (d) D = 0.5; (b) and (e) D = 2; (c) and (f) D = 20. In all cases e = 0.5, a = 1, 5 = 0.3, 
z = 1.5. 

FIG. 9. Time series generated by the noise-driven dynamical system (20). (a) D = 2, (b) 
D = 10. In all cases e = 0.5, a = 1, 5 = 0.1, z = 1.1. Note that in case (b) the control is so good 
that the output signal is almost periodic and is confined to a very narrow interval of values. Note 
also that in this case the noise intensity is 5 times larger! 
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FIG. 10. The explicit function (|23| ) can behave as a stochastic resonance system. SNR vs noise 
intensity {D) is shown for B = 24 and B = 1. Note that for B = 1 there is no maximum in this 
plot. 



FIG. 11. Chaotic time series generated by the continuous function (p4|). 
FIG. 12. SNR vs noise intensity (D) for the dynamics of system (16). Here the "noise" r](t) is 



defined as rj(t) = In [/(*)/ (1 - /(*))], where f(t) is given by Eq. (f24|). 

FIG. 13. SNR vs noise intensity (D) for the dynamics of system (p6|). The SNR is calculated 
from the time series generated by the function <j)(x = 0,t). 

FIG. 14. Time series produced by explicit functions that describe different stochastic resonance 
systems, (a) function (34), (b) system (35), (c) function (36), (d) function (39). 

FIG. 15. Typical time series for a bistable stochastic resonance system such as that described 
byEq. (16). 

FIG. 16. First-return map produced by function (39), which describes a stochastic resonance 
system with intrinsic noise, (a) Situation of stochastic resonance D = 0.16. (b) Situation out of 
stochastic resonance D = 0.3. 

FIG. 17. First-return map produced by the dynamics of variable Z n in the dynamical system 
(40)-(42). (a) Function g(x) possesses 1 local extremum. (b) Function g(x) possesses 100 local 
extremas. 

FIG. 18. I— V characteristic curves of two nonlinear circuits: (a) With some apropriately chosen 
input signal, this circuit can produce a very complex output signal (see the discussion in the text), 
(b) If the input signal is composed of a periodic signal and noise, the output signal will be less 
complex than the input signal. In fact, this will be a stochastic resonance system. 
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